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Abstract 

We study the transverse single spin asymmetry in Drell-Yan production in pA collisions with in¬ 
coming protons being transversely polarized. We carry out the calculation using a newly developed 
hybrid approach. The polarized cross section computed in the hybrid approach is consistent with 
that obtained from the usual TMD factorization at low transverse momentum as expected, whereas 
at high transverse momentum, color entanglement effect is found to play a role in contributing to 
the spin asymmetry of Drell-Yan production, though it is a 1/N^ suppressed effect. 


PACS numbers: 


I. INTRODUCTION 


Proton-nucleus (or deuteron-nucleus) collisions at RHIC provide an unique opportunity 
to study saturation/Color Class Condensate(CGC) physics. Many relevant observables in 
pA collisions in the forward rapidities region are excellent probes for accessing the saturated 
small X gluon distributions inside nucleus. The remarkable theoretical and experimental 
progress made recently in the held mostly focus on the spin independent observables, among 


ladron correlation at forward rapidities play 
Meanwhile, polarized proton-proton col- 


which the single hadron suppression and the di- 
important roles in studying saturation physics 
lisions at RHIC have made a big impact on the investigation of the nucleon’s spin structure. 
In particular, transverse single smn asymmetries(SSAs) phenomena in polarized pp collisions 
have ga^ned a lot of attenf.one g aa the afudy of SSAa aot only could help ua to n.ap out 
the three dimensional image of nucleon p|, but also greatly deepened our understanding of 
QCD and its associated factorization properties. 

On the other hand, scattering a polarized probe on a dense background gluon held in¬ 
side a large nucleus may provide a promising way of studyii^ the interplay of saturation 
ehects and transverse spin phenomena. The authors of paper have proposed to probe the 
saturation scale of nucleus by measuring SSAs normalized by that in pp scattering at low 
transverse momentum. It is also important to measure the SSA for prompt photon produc¬ 
tion in polarized pA collisions in order to distinguish diherent mechanisms for generating the 
SSA 5-7|. Furthermore, it has been shown that polarized observables are sensitive to the 

n 

slope of small x gluon transverse momentum dependent (TMD) distributions in k± space [8|- 
[i^ . Measuring transverse momentum dependence of SSAs thus may provide complementary 
information on small x gluon TMDs. 

Polarized pA collisions also present an advantage over unpolarized pA collisions and po¬ 


larized pp collisions in addressing one novel aspect of QCD: color entanglement ehect 




To generate the imaginary phase necessary for the non-vanishing SSAs, one additional gluon 
must be exchanged between the partonic hard part and the proton remnant. The interac¬ 
tions of this additional gluon and the valence quark from proton with the saturated gluon 
held inside nucleus lead to a very complicate color how structure which could give rise to 
color entanglement. Such ehect is the consequence of nontrivial interplay among the T-odd 
ehect, the coherent multiple gluon re-scattering, and the non-Abelian feature of QCD. In- 











vestigating SSAs in polarized pA collisions may shed new light on the study of generalized 
TMD factorization breaking effect that is caused by color entanglement. A polarized pA 
collisions program at RHIC is therefore extremely welcome 121. 

In this paper, we study the SSA in Drell-Yan lepton pair production at forward rapidities 
in polarized pA collisions. Due to the absence of hnal state interactions and fragmentation 
effects, the SSA in Drell-Yan process offers a very clean probe for the Sivers effect. The 
contribution from the Sivers effect to the SSAs has been well formulated in the context of 
TMD factorization and the collinear twist-3 approach. For the polarized pA collisions case, 
to incorporate the saturation effect, we carry out the calculation in a hybrid approach in 
which the nucleus is treated in the CGC framework while the collinear twist-3 formalism is 
applied on the proton side. Such hybrid approach has been recently developed to study the 
SSAs in prompt photon production and photon-jet production in polarized pA collisions 


Q. 


We notice that the SSA in Drell-Yan production has also been studied using a different 
hybrid approach ISj. As shown below, two different hybrid approaches yield the same result 
for this observable at low transverse momentum. 

In a more general context, the present work is part of the effort to address the inter¬ 
play between spin physics and saturation physics. Apart from the studies mentioned above, 
recent work in this very active held includes the study of the quark Boer-Mulders distribu¬ 


tion and the linearly polarized gluon distribution inside a large nucleus 


14i\ . The small x 


evolution equations for the linearly polarized gluon distributions were derived in Ref. 


The hrst numerical study of the linearly polarized gluon distribution was presented in 16 |. 


15|. 


Furthermore, the asymptotic behavior of transverse single spin asymmetries at small x was 


discussed in Ref. 


17 


18| . It has been shown that SSAs at small x are generated by the spin 


dependent odderon exchange whose size is determined by the anomalous magnetic moment 
of proton 1^. The quark Sivers function was computed in the Glauber-Mueller/McLerran- 
Venugopalan(MV) models 19|. The spin asymmetries in pA collisions have been investigated 
by going beyond the Eikonal approximation within the CGC framework 20|. 

The paper is structured as follows. In section II, we derive the spin dependent amplitude 
using the hybrid approach, including both soft gluon pole and hard gluon pole contributions. 
In section III, we present expressions for the polarized cross section in different kinematic 
limits and compare our results with that obtained from different approaches which are 
applicable in the corresponding kinematic regions. The paper is summarized in section IV. 












II. THE DERIVATION OF THE SPIN DEPENDENT AMPLITUDE 


In this section, we derive the spin dependent amplitnde for Drell-Yan production using 
the newly developed hybrid approach. We start by briefly reviewing the CGC calculation 
for unpolarized Drell-Yan cross section in pA collisions. 

The dominant production mechanism for Drell-Yan virtual photons at forward rapidities 
is Compton scattering qg —)■ 7 *g. We £x the relevant kinematical variables and assign 
4-momenta to the particles according to 


q{xP) + g{x'P + ki_) —^ 7*(V) + 


( 1 ) 


where P^ = P~n^ and P^ = P~^pi^ with and being the commonly dehned light 
cone vectors, normalized according to p ■ n = 1. The Mandelstam variables are dehned as: 
S = {P + Py, T = {P — IqY and U = (P — l^*)"^. The invariant mass of the produced lepton 
pair is denoted as It is worthy to mention that x'^P + k± is the total momentum 

transfer via multiple gluon re-scattering. 

The calculation for Drell-Yan virtual photon production in unpolarized pA collisions is 
rather similar to that for prompt photon production, and has been done within the CGC 
framework a decade ago 2l(]. The key ingredient of this calculation is resumming multiple 
gluon re-scattering into Wilson line which is a path-ordered gauge factor along the straight 
line that extends in x'^ from minus inhnity to plus inhnity. More precisely, for a quark with 
incoming momentum I and outgoing momentum I + A:, the path-ordered gauge factor reads. 


27r6{k+)nf^[U{k^) - {27rf 6{kP] , 


( 2 ) 


with 


U{k^) 




and 

U{x±) = VP9p^dx+A-ix+, x^yt ^ 


( 3 ) 

( 4 ) 


where t is the generators in the fundamental representation. With this calculation recipe, it 
is straightforward to obtain the cross section for unpolarized Drell-Yan lepton pair produc¬ 


tion 


2l|, 


Pa 


dM‘^(PL*±dy 


q 




dx 


X 


d'^k± H{k±,p*_i_,z)xfq{x)x' GDp{x'k±) (5) 






where y is the rapidity of the virtual photon. In the above formula, fq{x) is the integrated 
unpolarized quark distribution from proton, and x'gGDp{x'g, kj_) is the dipole type gluon 
TMD, defined as. 


xfinp(x,M) - ^ 


it,.1 

(27r)2 N, 




( 6 ) 


The hard part H{k^, z) is given by. 






[{lj*± — zk±Y + e\j\ \l. 
1 1 


+ 


n 2 


1'^* 


+ e 


M 


{Y*± — zk±y + e 


M 


( 7 ) 


where z = Y* ■ n/{xP ■ n) is the fraction of the incoming quark momentum xP carried by 
the virtual photon, and e\^ = {1 — z)!^"^. Y*± is the virtual photon transverse momentum. 
The connections between the above Drell-Yan cross section derived in the CGC framework 
and those from the TMD factorization and the collinear factorization have been discussed 


m paper 


13 |- 


We now move on to derive the spin dependent amplitude of Drell-Yan production. To 
generate the spin asymmetry, one additional gluon must be exchanged between the active 
partons and the remnant part of the polarized proton projectile. In the collinear twist-3 


approach, t 
function (22 


re associated soft part is described by the three parton correlator: the ETQS 


TF,q{xi,x) = 


dyi dy. 


47r 


2 +i{x-xi)P+y 


x{P,SM,(Oh*gpx,''<’F^+(y-)i,,(yp\P,Sp) 


( 8 ) 


where we have suppressed Wilson lines. S± denotes the proton transverse spin vector. 

As mentioned above, we derive the spin dependent amplitude using a hybrid approach 
which is formulated in the covariant gauge. The additional exchanged gluon is longitudinally 
polarized in such covariant gauge calculation. Unlike a quark scattering off the classical 
color held of nucleus, the multiple scattering of this longitudinally polarized gluon with 
the background gluon held of nucleus can not be simply described by a Wilson line in the 
CGC formalism. Instead, the expression for the gauge held created through the fusion of 
the incoming longitudinally polarized gluon from the proton and small x gluons from the 















FIG. 1: Diagrams contributing to the spin dependent Drell-Yan lepton pair production amplitude. 
Different symbols indicate different parts of the classical gluon field. A black dot denotes A^eg or 
Aging, while a cross surrounded with a circle denotes Aa- Diagrams Fig.la, Fig.lb, Fig.ld, Fig.le, 
Fig.If, and Fig.li-Fig.lp generate the soft gluon pole contribution. Diagrams Fig.la-Fig.ld and 
Fig.li-Fig.ll give rise to the hard gluon pole contribution. 


nucleus takes a quite couiplicate form, and contains both singular terms (proportional to 
and regular terms |24| . 


A^{q) = A^(q) + 5^-A-,(q) . 


The regular terms are given by 


/IM _ /IM 
■^reg 




(PPA 


,P_l) U{ki_) - 


q^ + iq+e J (27r)2 

+ regig) “ (2vr)^5(fcx) 


I PpjPl 

J p\ 


(9) 


( 10 ) 


where Ppipi) is the color source distribution inside a proton, and A^ is the gauge field created 
by the proton alone. In the second term of the formula [101 is the momentum carried 































by the incoming gluon from the proton and k± dehned as k± = q± — p± is the momentum 
coming from the nucleus. For the polarized case, there exists a correlation between the 
transverse momentum p± and the transverse proton spin vector S±. Such a correlation is 
described by the ETQS function, and leads to a SSA for direct photon production jo] as well 
as Drell-Yan lepton pair production. The four vectors C^{q,p±) and are given by the 

following relations, 

1,2 

C'f^(g,P±) = 

y T oc 

p. 




+ le 


C\j{(i,P2) = -2pj 




T 


( 11 ) 

( 12 ) 


+ le 

where the subscript 'reg' indicates that the corresponding term of does not contain any 
6{x^) when expressed in coordinate space. Here, we specihed the q~^ pole structure according 
to the fact that this term arises from an initial state interaction. It is crucial to keep the 
imaginary part of this pole in order to generate the non-vanishing spin asymmetry. The 
notation pj_ is used to denote four dimension vector with p^ = —p^^. U{k±) and V{k±) are 
the Fourier transform of Wilson lines in the adjoint representation, 

__ r r /•-i-oo 


V{ki_) = / (fx±e 


ik±-x± 


Vexp 


^9 


■ 9 


’ —OO 

r+oo 


dz~^A^{z~^, x±) ■ T 
dz^A'^{z^, xj,) ■ T 


(13) 

(14) 


where the T are the generators of the adjoint representation. The singular terms reads. 


AingiQ) = 




_ f d'^p^ 

q+ + ie J {2Try 


V{k^) - {27iY6{k^) 


PpjPi 

p\ 


(15) 


The peculiar Wilson line V differs from the normal one f/ by a factor 1/2 in the exponent. It 
appears to be a generic feature that all terms containing V cancel eventually when computing 
a physical observable jb, 10, 24, 251. 

Following the method outlined in Ref. 25|, one has to calculate the contributions from 
the regular terms and the singular terms separately. In the prompt photon production case, 
the imaginary phase necessary for non-vanishing SSAs is generated from the soft gluon pole 
while for Drell-Yan production, the imaginary phase arises from both the soft gluon pole 
and the hard gluon pole due to the existence of an additional hard scale M. 

The derivation of the spin dependent amplitude which contains the soft gluon pole con¬ 
tribution is very similar to that presented in Ref. jb]. The hnal expression for this amplitude 























takes form, 


Msgp = 


d^P± PpAp^) f dki d^kii_ 


(27r)2 pI 


(27r) 


U{lq 


X <; ^Sf{x,P + + ki)iU{k,^) 

yq — rCi) + *6 

(Pu{q - p ki)iU{ki±)SF{xiP - 


(g - A;i)2 + ie 

(puiq - kup±). 


rSpixiP — /p + ki)^SF{xiP + ki)7^U {ki±) 


+ ^eg 


iq-k,)^ + ie'' 
xu{xiP) U{k± — ki±) — {271)"^6{k± — ki±) 
d‘^P± Pp,a{P- 


ha 


(27r)2 pI 


d x±e 


ik±-x± 


Xu{lq 


4Sf{xiP - /p)^ + ^Sf{xiP + q)i f, 


XgP + ie 


fU{xF)u{xiP) \ U{xf) - 1 (16) 


ha 


The hard gluon pole is generated when the quark propagator Sf{xiP—1^*) = 


+ie 


goes on shell. It provides a phase proportional to S{xg — Xg) where Xg = x — xi = Xg = 
X — M‘^/2P ■ Diagrams Fig.la-Fig.ld and Fig.li-Fig.ll give rise to the hard gluon pole 
contribution. The corresponding spin dependent amplitude reads, 


Mhgp = -ieg'" 


d'^Px PpAP-l) f dk^ (Pki± ^ 

U[lq) 


(2^)2 pI 


(2ir) 


X 


I (Puiq - V + k,)iU{k,F)SF{x,P - 

[ [q- ki)^ + te 


(puiq - ki,p±). 


P-PSf{xiP — + ki)^SF{xiP + ki)7^U {ki±) 


+ ieg 


{q-kiY + ie^ 
xu{xiP) U{k± — ki±) — {27r)‘^6{k± — ki±) 

d^px Pp,a{p. 


ba 


(27r)2 p\ 


d x±e 


ik±-x± 


X <j u{lq){x±)u{xiP) \u{x±) - 1 
' XqF + le L . 


ba 


+iu{lq)p/PS f{xiP — + k)'tji \U (xj,) — 1] Sf{xiP — l^*)^u{xiP) 

+iu{lq)7^ [U{x±) — 1] — k)plPSF{XiP — l^*)^u{XiP)^ (17) 


The last two terms in the above formula come from the Ap part of Areg in Fig. lb and 
Fig.lc. Note that Fig.lc contains no soft gluon pole, while the soft gluon pole contribution 
from Fig.lb cancel out with its conjugate diagram. With these two derived amplitudes, it 
is difficult to compute the full twist-3 polarized cross section. However, in the kinematic 



regions where the collinear approach or the TMD factorization approach is applicable, the 
calcnlation can be greatly simplihed, snch that we can make comparisons among different 
formalisms. 


III. THE SSA AT LOW AND HIGH TRANSVERSE MOMENTUM 


The SSA of Drell-Yan prodnction can be described in the context of the collinear higher 
twist factorization approach at high transverse momentnm and the TMD 

factorization framework at low transverse momentnm In this section, we com¬ 

pare onr hybrid approach with these two formalisms in the corresponding kinematic regions. 

We hrst extrapolate the fnll amplitndes to high transverse momentnm region by Taylor 
expanding the hard parts in terms of ki±. Repeating the power connting analysis in Ref. jb] 
and neglecting power snppressed contribntions, the fnll amplitndes can be dramatically 
simplihed to. 


M 


SGP 


M 


HGP 


-leg 


(Pp± Pp,a{P-^ 

(27r)2 pI 


(Px±e 


ik^-x±_ 


U{lq) 


X 


(pL{.q,PA.)SF{,XiP - + ^Sf{XiP + q)(pL{.q,Pp) .b 


-leg 


(Pp± Pp,a{Pl 
{2fY pI 


-|- ie 
/ (fx±e 


rU{x±)u{xiP) U{x±) 


ba 




X i u{lq) 


<pL{q,P±)SF{xiP - lr)Lb 


-\- ie 


fU{xF)u{xiP) U{xf)-1 


ba 


-iu{lq)jlPS f{xiP — Iry* k)r^ — 1 ] Sf{xiP — l^*)^u{xiP) 

-iu{lq)l^ \U{Xf) — 1 ] — k)q/PSF{XiP — Iy)^u{XiP)^ 


(18) 


where (pL is the well known effective Lipatov vertex for the prodnction of a glnon via the 
fnsion of two glnons, and given by. 


(pL _ (pu _ 

q^ + ie q^ + ie g+ -|- ie 


(19) 


The next step is to fnrther expand the hard part in terms of k^ and keep the leading power 
contribntion following the method ontlined in Section 3.3 in Ref. jb]. After having done 
so, it becomes evident that the hard coefficient calcnlated from the above amplitnde is the 


same as that compnted in the standard collinear twist-3 factorization 


2b I. On the other 


hand, nsing the Fierz identity, the soft part from the nnclens side, i.e. Wilson lines, can 
















be reorganized and expressed into two parts: unpolarized gluon distribution and the novel 
gluon distribution G 4 . Collecting the hard gluon pole contribution and the soft gluon pole 
contribution, we eventually obtain the following polarized differential cross section, 


dAa 


AirOig^m cij3 


dM‘^d?U^dy 3N^SM^ 47r2 




X 


<? 
d 


dxdx'^^ . ^ 1 

- -dis + t + u — M ^)—^ 

X x' —u 


HBorn ( TF,q{x,x) - X—TF,g{x,x) ) + NsTF,g{x,x) 


H' 


HGP 


+ TF,g{x -Xg,x) 

where the hard coefficients are given by, 

HBorn 

iV, = 


[Gnp{x'g) - G 4 K)] 
[Gopix'g) - GA{x'g)\ + H^FGpGDp{x'g) 


( 20 ) 


iV2 

S 

i 2 M‘^u 

IVI-l 

i 

s si 




[M^ - 2MH + 


H 


iV 2 - 1 sP 

(M2 - if + 


HGP 


m 


si^ 


s (M 2 - if + 

S + U sf 2 


H 


HGP 


( 21 ) 

( 22 ) 

(23) 

(24) 


with s, u, t being the normal partonic Mandelstam variables. In the collinear limit lA, S> kf 


they can be expressed as. 


s = 




u = 


p, 

S* 


t = 


P, 


(25) 


{l-z){l-z) 1-z “ 1-z 

where 5 is the longitudinal momentum fraction of the incoming gluons x'gP carried by the vir¬ 
tual photon. In the above formula, x'gGFip{x'g) and x'gGfxf kj_) are the integrated gluon dis¬ 
tributions dehned as x'gGrtpix'g) = / d‘^kj_x'gGFip{x'g, k^), x'gGfx'g) = / dPkFx'gGfxfkf), 
respectively. The gluon distribution Gfxf kf) possesses an unique Wilson line structure, 

%Jy2.e“-(--!'ji(Tr4C7(.x)lTr,[C7l(!,x)l>.; (26) 


It has been shown in the MV model [27| that G 4 is suppressed by the power of 1/A^"2 as 
compared to the unpolarized dipole distribution jsl. If we neglect these l/iV2 suppressed 
terms that essentially arises from color entanglement effect, it is easy to see that the polarized 
cross section computed in the standard collinear twist-3 approach can be recovered from our 
hybrid approach. 
























To compare with the result from TMD factorization, we have to further extrapolate the 
above result to the moderate transverse momentum region, << << where TMD 

factorization is applicable. This can be easily done by expanding the delta function. 


6{s + i + it — M^) = 


(5(1 — z) (5(1 — 5) , , , 

- - + 7 T-+ <5(1 - z)5{l - z)\nj^ 


(1 


(1 


7 *± 


(27) 


where only the hrst term proportional to (5(1 —; 2 ) gives rise to the leading power contribution. 
In the limit z ^ 1, one has |f| << !«! ~ s, which implies Xg —?■ 0. The fact that the hard 
gluon pole degenerates with the soft gluon pole in this kinematic region allows us to combine 
two contributions together. Substituting the above expansion into Eq. [20l one ends up with. 


dAa 






dM‘^d?l^*±dy 


3N^SM^ 47r2 


/4 


X 


-^(5(1 - z)TF,q{x, x)GDpix'g) [z^ + (1 - zf] 


X X' 


(28) 


which agrees with Eq.31 in Ref. 2^. As expected, all terms arising from the color entangle¬ 
ment effect cancel out in the kinematic limit we consider. Therefore, our hybrid approach 
is consistent with TMD factorization at moderate transverse momentum. 


However, TMD factorization applies in a broader kinematic region: -C M^. On 

the other hand, the hybrid approach is valid as long as A^cd ^ To demonstrate the 
complete equivalence between the two formalisms in the overlap region "C -C M^, 

it is necessary to keep k± hnite when computing the hard coefficients. This makes the 
evaluation of the polarized cross section much more involved. Nevertheless, through the 
explicit calculation, we verify that the SSA does not receive the contribution from the initial 
state interaction due to the complete cancelation between the soft gluon pole contribution 
and the hard gluon pole contribution in the low transverse momentum region. We are 
thus left with the hard gluon pole contribution from diagrams Fig.lb and Fig.lc. The 
corresponding amplitudes are give by the last two terms in Eq. [T71 We further found that 
the hnal state interaction shown in Fig. lb only yields the power suppressed contribution at 
low transverse momentum. The only remaining piece is the hard gluon pole contribution 
from diagram Fig.lc which can be computed following the standard procedure. The fact 
that the entire surviving contribution is just due to hard pole at low transverse momentum 
has also been observed in Ref. 29|]. At this step, we would like to mention that the soft 


part associated with the diagram Fig.lc only contains the regular Wilson line structure. In 












the end, the polarized cross section at low transverse momentnm can be nicely cast into the 
following compact form, 


dAa 


Airalm as 1 


dM‘^d%*^dy 3N^SM^ 2 


X 




X' 


kl 




(29) 


with A(A:x, ^ 7 *±, 5) being dehned as. 




V-lIV-L “ ^-lI 


— /cj 


{I - z)l‘i^^ +k±Y \l^*±-kx\ 


(30) 


The transverse momentnm carried by small x glnons is of the order of the saturation scale 
Qs- At moderate transverse momentum S> l'y*± 3> Qs ~ k±, one can Taylor expand 
the hard coefficient dA{k_L,l-y*±, z)/dl^,j_ in terms of kj_. By keeping the nontrivial leading 
term and neglecting the terms suppressed by the power of k_i_/lY±, Eq. |28] can be readily 
recovered from Eq. [211 

As mentioned above, the SSA for the Drell-Yan production at low transverse momentum 
also can be described in the TMD factorization approach. The corresponding polarized cross 
section reads. 


dAa 


Anal 


dM‘^d%*^dy 3N^SM^ 


X 


E 




PX-p'A_)ftT,q{x,Pl.)fq{x',p'^) (31) 


where and fq{x',p'j_) are the quark Sivers function and the unpolarized anti-quark 

distribution from the target nucleus, respectively. At small x, the anti-quark distribution is 
dynamically generated through gluon splitting process j^ . 

n r dx' r d'^k , 

f,i^'^P'±) = ^ 2 ] ^GMx'^,k^)A{k^,p'^,z) (32) 

As argued above, the typical small x anti-quark transverse momentum p'j_ = — pj_ is 

of the order of Qs and much larger than the incoming quark transverse momentum p± ~ 
Aqcd- After substituting Eq. [32] into Eq. [31] and carrying out the integration over p'j_, we 
make Taylor expansion for A{kj_, Q*± — p±, z) in terms of p± and keep the linear term. To 
make the connection between Eq. [H] and Eq. [21] one further use the well known relation 















Tp^q = / d‘^p±^f^^q{x,p±) for the Drell-Yan process [30|. It is then straightforward to 
reproduce Eq. |29] from Eq. [3ll 

We thus conhrm that the hybrid approach agrees with TMD factorization for the Drell- 
Yan production in the full overlap kinematical region where both the formalisms are appli¬ 
cable. This agreement provides strong evidence that the hybrid approach is complete and 
self-consistent. We notice that the equivalence between the small x formalism and the TMD 
factorization approach in describing the same observable has also been verihed in Ref. Q- 


IV. CONCLUSION 


Color entanglement effect is usually believed to be absent in the Drell-Yan process in the 
context of TMD factorization because of the simple color flow structure, though a complete 
consensus has not yet been reached [Ml]- With the help of newly developed hybrid approach, 
we are able for the hrst time to provide a non-trivial check on this statement for the SSA 
case. To be more precise, we take into account one additional gluon exchange from polarized 
proton and resum gluon scattering to all orders on nucleus side using the hybrid approach. 
It has been shown that in the polarized cross section, all terms arising from the color en¬ 
tanglement effect drop out at low transverse momentum due to the systematical cancelation 
between the soft gluon pole contribution and the hard gluon pole contribution. However, 
at high transverse momentum, the polarized cross section computed in the hybrid approach 
differs from that obtained from the standard collinear twist-3 approach by some additional 
contributions whose emergence can be attributed to color entanglement. Such novel color 
entanglement effect in principle could be studied at RHIC [12|, though it is found to be 
l/N^ suppressed and thus very small. 
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